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ABSTRACT 


PROJECT  BRI3HAM  was  a Department  of  Defease  data 
collection  effort  wherein  890  kHz  wide  samples  of  the 
HF  spectrum  were  acquired  using  a 25  Hz  sampling  rate 
with  2.3  kHz  resolution  for  a period  of  up  to  2.4 
minutes.  The  main  tacust  of  this  paper  has  oeen  rhe 
visual  examination  of  the  data  on  a graphics  display 
with  the  idea  of  identifying  drastic  and  unexpected 
changes  in  either  individual  signals  or  tne  entire 
wavefront;  those  changes  being  due,  at  least  in  part, 
to  signal  transformation  resulting  from  various 
ionospheric  phenomena.  In  addition,  an  in-depth 
statistical  study  was  conducted  on  several  of  the 
specific  signals  witn  the  hope  of  aiding  in  the 
formulation  of  an  algorithm  to  efficiently  automate 
the  receiving  system's  sensitivity  settings. 
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I. 


OBJECTIVES 


The  objectives  of  the  project  were:  (1)  to  continue  the 
visual,  or  graphical,  analysis  of  the  BRISHA'l  data  (the 
initial  effort  of  which  was  reported  in  Ref.  1)  in  an 
attempt  to  verify  anomalous  signal  behavior  already 
documented  and  to  discover  any  new  anomalies  wnica  might  ne 
present;  (2)  toward  that  end,  to  malce  whatever  caanges  were 
necessary  in  both  tne  programs  and  data  handling  procedures 
to  meet  either  new  system  constraints  or  to  enhance  over-all 
project  efficiency;  and  (3t  to  analyze  the  data  from  a 
statistical  standpoint  (via  two  entirely  new  computer 
programs)  in  an  effort  to  determine  how  often  and  by  what 
amount  (in  decibels)  several  different  classes  of  signals 
change. 

The  format  of  the  data,  particularly  that  of  the  2. '4 
minute  time  duration  limitation  (although  there  were  several 
time-wise  contiguous  data  sets) , made  long  term  analysis, 
botn  visual  and  statistical,  impossible.  In  addition,  the 
daily  intercept  times  for  the  data  sets  precluded  any 
significant  effort  being  made  to  correlate  signal 
cnaracterist ics  from  even  approximately  the  same  time  of  the 
day. 


II. 


INTRODUCTION 


It  is  well  documental  that  as  radio  signals  traverse  the 
ionosphere  (including  signals  bounced  off  of  one  of  its 
layers)  they  ara  often  cnanged  in  a manner  that  goes  beyond 
simple  attenuation  due  to  the  distance  travelled.  Both  the 
make-up  of  the  ionosphere  itself  (how  many  ions  are  present 
at  any  one  time,  tna  alusive  Sporadic-E  phenomenon, 
niynt-day  asymmetry  and  the  geomagnetic  non-reciprocity  to 
name  but  a few)  and  the  more  dramatic  disruptions  of 
ionospheric  behavior  (solar  flare  caused  magnetic  storm 
activity,  polar  cap  absorption  and  tae  auroral  displays  are 
three)  team  with  the  more  common,  every  day  pnenomena  of 
"Faraday  Rotation"  and  noise-producing  lightning  storms  to 
alter  any  signal  with  which  any  or  all  of  taem  come  in 
contact.  The  results  can  ds  severe  attenuation,  pronounced 
amplification,  "ducting",  "whistlers",  and  fading  (sometimes 
called  the  Dellinger  Sffect) . References  2-5  were 
excellent  sources  for  ionospheric  phenomena  information. 

It  was  felt  that  the  application  of  graphical  techniques 
to  signals  which  had  traversed  the  ionosphere  aignt  add  some 
insight  into  the  problem  of  exactly  what  happens  to  tnese 
signals. 

For  the  most  part,  the  visual  presentation  techniques, 
shown  in  various  stages  in  Fig.  1,  us=d  in  this  study  were 
the  same  as  those  presented  in  Re£^  1.  However,  some 
changes  have  been  made  and  will  be "discussed  in  detail  in 
the  appropriate  section.  as  in  Ref.  1,  interesting  sets  of 
scans  have  been  reproduced  as  graph  plots  (using  the 
VERSATEC  plotter)  for  further  analysis  and  presentation  in 
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this  paper.  All  thirty- fair  of  the  previously  unanalyzed 
complete  lata  sets  were  studied  with  varying  degrees  of 
success  concerning  anomalies  (or,  least,  perceived  to  be  soi 
found  therin. 

The  swept-tuned  receiver  and  associated  equipment  used 
to  acguire  and  record  the  data  sets  is  an  operational 
receiver  eguipment  used  for  otner  than  experimental  data 
gathering.  Its  proper  operation  requires  frequent  tuning  of 
numerous  sensitivity  adjustments.  Ho*  often  and  by  what 
amount  these  attenuators  and  thresholds  are  cnaaged  have  a 
profound  effect  on  the  performance  of  the  eguipment. 
Therefore,  programs  MIN/dAX  1 and  MIN/MAX  2 were  written  in 
an  effort  to  accumulate  enough  statistical  data  from  the 
sets  to  warrant  the  recommendation  of  some  appropriate 
algorithm  for  controlling  the  attenuator  adjustments. 
Section  IV  of  this  paper  will  address  itself  to  chat  topic. 


A.  NATURE  OF  THE  DATA 


The  data  was  gathered  using  a swept-tuned  receiver  with 
either  an  omnidirectional  or  12  degree  beam  antenna.  The 
data  sets  were  produced  oy  sampling  the  receiver  output  at  a 
25  Hz  sampling  rate  as  ic  scanned  downward  through  an  390 
kHz  band.  Each  data  set  consists  of,  at  most,  3597  scans 
(each  scan  being  composed  of  318  contiguous  2.3  kHz  wide 
bins  moving  down  in  frequency  as  the  bin  numoer  increases) 
placed  on  7194  cards.  Tied  one  gives  amplitude  information, 
in  2 dB  high  quantization  levels,  for  bins  1 througn  150 
while  information  on  bins  161  through  318  is  found  on  card 
two.  The  quantization  levels  are  given  in  5-oit  binary 
words  with  even  numbered  bin  information  found  in  rows  11 
through  3 while  odd  nunmered  bin  information  is  in  rows  5 
through  9.  In  addition,  each  card  a as  the  card  numoer  and 
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set  number  punched  in  row  12  for  accounting  purposes.  With 
each  set  there  is  a paraaeter  card  which  delineates  the  date 
and  rime  of  signal  acquisition,  base  frequency,  antenna 
configuration  and  attenuation  applied. 

Of  approximately  150  data  sets  recorded,  forty-four  were 
made  available  for  this  study.  Figure  2 is  a table  of 
pertinent  information  concerning  the  forty-two  complete  sets 
analyzed  in  Ref.  1 and  this  paper. 

B.  CHANGES  IN  THE  GRAPHICAL  ANALYSIS  PROCEDURE 


Equipment  reconfigurations  since  the  graphical  analysis 
procedure  was  formulated  and  reported  in  Ref.  1 and  efforts 
to  make  the  data  .handling  more  efficient  demanded  several 
changes  both  in  procedure  and  in  the  two  main  analysis 
programs.  Anomaly  A and  Anomaly  B.  Specific  information 
concerning  the  wordings  of  the  two  programs  may  oe  found  in 
Ref.  1,  particularly  in  section  II  B. 

First  of  all,  the  sheer  number  of  binary  cards  to  be 
transcrioed  onto  magnetic  tape  (nearly  a quarter  of  a 
million)  obviated  the  card-to-tape  technique  used 
previously.  A much  faster  and  far  more  reliaoLe  procedure 
was  initiated  wherein  tie  I3M  363/57  system,  belonging  to 
the  W.  R.  Church  Computer  Center,  was  used  after  normal 
working  nours.  The  cards  ware  first  read  onto  9-track  tapes 
(800  BPI,  even  parity)  and  then  transcrioed  to  the  7-track 
tapes  (556  BPI,  odd  parity)  required  for  the  XDS-9300 
computer.  Partial  dumps  were  made  of  all  sets  from  both  the 
9-  and  7-track  tapes  to  ensure  that  the  cards  read  onto  the 
tapes  properly  and  that  tie  tape-to-tape  translation  ns 
successful.  The  programs  used  to  read  the  cards  onto  tapes, 
translate  the  9-  to  7-tract  and  dump  tne  tapes  for 
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quality  control  cm  be  found  at  the  end  of  this  paper. 

This  nev  system  of  data  loading  resulting  in  the 
requirement  to  cnange  the  sections  of  the  Anomaly  A and 
Anomaly  3 programs  refering  to  data  read-in.  This  was 
accomplished  with  little  difficulty. 

Next,  a metasymbol  subprogram  was  added  to  each  of  tn? 
analysis  programs  to  enaala  entire  files  (data  sets)  to  be 
skipped  when  a specific  file  required  was  located  other  than 
first  on  a tape.  Ey  calling  the  routine  "P0R5CN"  a 
keyboard-entered  number  of  files  (NFILE  = ) would  be  passed 
over.  This  proved  to  be  extremely  valuaole  as  a time  and 
labor  saving  device.  The  coutine  ''BAKSCN'1  was  also  availanle 
for  Dacking  over  entice  files  but  cnronic  tape  drive 
malfunctions  negated  the  use  of  this  even  more  valuaole 
tool. 

One  final  major  program  change  was  dictated  by  a 
computer  center  equipment  reconfiguration  wherein  tne 
CALCOMP  563  plotter  was  no  longer  available.  The  new  system 
to  be  used  was  the  PDP-11  computer  in  conjunction  with  a 
V ERSkTEC  plotter.  This  cnange  required  modifications  of  tne 
analysis  program  sections  dealing  witn  tne  processing  of  the 
output  data  (desired  plots  from  the  AGT-10  screen),  found  at 
the  very  end  of  tne  main  program,  and  the  generation  of  a 
new  suoroutine  , called  the  Header  Subroutine,  which  packs 
the  header  data  into  the  proper  format  for  the  plotting 
unit.  In  addition,  a PDP-11  plot  package  (used  to  interface 
the  data  tape  wit-h  the  plotter)  had  to  be  extensively 
modified  to  meet  the  demands  of  this  anomaly  data.  Listings 
of  the  revised  Anomaly  & and  B programs,  the  metasymbol 
program  for  calling  "FOasCN"  and  "BAKSCN"  and  the  PDP-11 
plot  routine  (called  "aa.c"  for  Anomaly  A plots  - this 
routine  alone  is  listed  since  tne  one  utilized  for  the 
Anomaly  B graph  plots  differs  only  in  a few  of  the 
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approximately  433  lines  comprising  "aa.c")  can  ds  found  at 
tne  back  of  this  paper.  Subroutines  3INP,  3INPUT  , FNS  and 
VCD  were  not  changed  and,  nence,  ace  not  listed  in  this 
paper  but  may  be  found  in  Ref.  1. 

C.  ANOMALY  A AND  B PLOTS 

The  plots  which  have  resulted  from  the  graphical 
analysis  portion  of  tai  s study  reguire  some  prior 
explanation.  Due  to  the  data  nandiing  restrictions  imposed 
by  the  XD3-9300  computer,  only  a small  portion  of  the  data 
in  each  set  could  be  viewed  at  any  one  time.  This  resulted 
in  some  difficulty  in  visualizing  what  could  be  long-term 
(time  wise)  or  wide  (freguency  wise)  phenomena  and  demanded 
that  each  set  be  analyzed  several  tines  before  the  entire 
set  had  been  viewed. 

The  Anomaly  A plots  consist  of  twenty  scans  (signal 
wavefronts  sampled  0.34  seconds  apart  - or  further  if  scans 
were  deliberately  skipped  to  view  long-term  effects)  with 
the  most  recent  scan  in  time  falling  at  the  bottom  of  tne 
picture.  The  scans  could  be  brought  onto  the  screen  in  any 
number  up  to,  and  including,  twenty  at  a time.  \s  shown  in 
Fig.  3,  the  vertical  axis  is  marked  off  in  5 da  increments 
which  apply  only  to  tne  bottom  scan  (called  "base  scan"). 
The  horizontal  axis  is  divided  into  ten  23  kHz  wide  sections 
each  composed  of  ten  freguency  bins  2.8  kHz  wide  for  a total 
of  100  freguency  bins  and  230  kHz.  Tais  width  can  be  varied 
to  other  than  100  bins  if  desired.  The  leftmost  oin  (termed 
the  "base  bin")  is  selected  by  keyuorad  input,  thus  allowing 
the  analyst  to  slide  up  and  down  the  data  set 
(freguency- wise)  at  will.  The  scale  along  the  horizontal 
axis,  then,  consists  of  increments  of  ten  bins  starting  with 
the  "base  bin"  number  at  the  origin.  The  term  "oase  freg" 


refers  to  the  base  freguency  of  the  data  set  and  decreases 
in  2.8  kHz  increments  as  the  bia  number  increases.  Shading 
has  been  provided  on  some  of  the  plots  to  highlight  the 
suggested  anomaly. 

The  Anomaly  B plots  are  a combination  of  the  Anomaly  A 
plot  technique  and  a modification  wnich  enables  tna  analyst 
to  view  selected  bins  (signals)  in  the  time  vs.  amplitude 
domain.  At  the  bottom  of  each  plot  are  five  scans  (Anomaly 
A style)  used  to  aid  in  the  selection  of  the  bins  to  be 
plotted  in  the  upper  traces  (see  Fig.  4).  in  addition  to 
the  restricted  number  of  scans  to  be  seen  at  one  time  (only 
5) , computer  limitations  required  tnat  the  number  of 
frequency  bins  presented  be  reduced  to  eignty.  Other  than 
those  changes,  discussion  of  the  Anomaly  A plots  is  also 
valid  for  the  lower  set  of  Anomaly  3 plots.  The  upper  set 
of  Anomaly  B plots  is  a collection  of  five  keyboard-selected 
signals  cnosen  from  among  the  scans  below.  The  traces  are 
numbered  one  to  five,  from  top  to  bottom  respectively.  The 
horizontal  axis  is  a time  line  composed  of  eighty  data 
points  at  least  3.04  seconds  apart  (depending  upon  the 
number  of  scnas  skipped,  if  any)  witn  tne  oldest  point  at 
the  extreme  right  side  of  the  trace)  with  the  hasi  marks  on 
the  vertical  axis  refering  to  13  d3  amplitude  levels. 

As  seen  to  some  extent  in  the  Anomaly  A plots  and,  more 
markedly,  in  the  Anomaly  B plots,  tne  scans  look  fractured 
in  appearance.  This  was  not  the  result  of  intentional 
programming,  rather  it  was  the  product  of  the  new 
PDP-1 1/VERSATEC  plot  routine  connination. 
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GRAPHICAL  ANALYSIS 


A.  DISCUSSION 

Prior  to  initiation  of  tne  grapnical  analysis  process 
via  programs  Anomaly  A and  3,  an  enumeration  of  many  of  the 
possible  "anomalies"  expected  was  made.  This  was  not  to 
imply  either  that  signals  which  had  Deen  markedly  affected 
by  the  ionosphere  would  only  exhibit  these  graphical 
manifestations  of  anomalous  performance  or  that  all  forms  of 
signal  anomalies  would  be  seen  oy  the  analyst.  Tais 
enumeration  procedure  was  merely  an  attempt  to  'visualize" 
what  effects  those  siga al-chaagiag  ionospheric  phenomena 
discussed  in  the  introduction  might  have  on  the  data  sets 
yet  to  be  analyzed.  Althouga  no  attempt  was  made  to  make 
this  study's  "anomalies"  conform  witn  those  identified  in 
Ref.  1,  there  has  been  an  attempt  to  verify  any  previously 
seen. 

The  graphical  effects  on  the  signal  waveforms  expected 
were: 

1.  a sudden  enhancenent/depression  of  one  or  more 
signals  (long  or  short  term); 

2.  an  abrupt  discontinuity  across  an  entire 
wavefront  or  poction  thereof; 

3.  a signal  "bum?"  or  "valley"  moving  eitaer  up 
or  down  in  freguency  as  the  wavefront 


1 9 


r ^ 


progresses  (Pig.  4,  Ref.  1 refers); 

4.  a shift  in  freguency  of  one  or  more  signals 
(either  short  term,  long  term  or  aperiodic  in 
nature)  ; 

5.  an  unexplainable  enhancement  of  one  or  more 
signals  with  a concomitant  depression  of  one 
or  more  signals  at  the  same  time  in  the  same 
data  set; 

6.  a complete  loss  of  continuity  of  the  signal 
wavefront  for  variable  periods  of  time; 

7.  unexpected  variations  in  one  or  more  signal 
contours  as  the  wavefront  progresses; 

8.  any  out- of-tha- ordinary  graphical  implication 
of  anomalous  performance. 

Effects  1,  2,  4,  5,  and  3 were  discerned  in  the  present 
study  while  several  of  these,  plus  effects  3,  5,  and  7,  were 
reported  in  Ref.  1. 

d . ANOMALY  A ANOMALIES 

As  shown  in  Ref.  1 and  verified  in  this  paper,  most  of 
the  data  sets  (whan  viewed  on  the  graphics  terminal) 
exhibited  a great  measure  of  uniformity  in  signal  makeup  in 
tne  time  domain  and  did  not  yield  any  obvious  graphical 
anomalies.  In  fact,  fully  fifty  per  cant  of  tne  BRIO. 8 AM 
data  sets  tailed  to  provide  any  reportaole  irregularities  in 
signal  pattern  with  tne  otner  fifty  per  cent  giving  rairly 
uniform  signal  wavefronts  for  most  of  the  time  recorded. 
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The  latter  group  did,  however,  exhibit  anomalous  (or,  at 
least,  assumed  to  be  so)  characteristics  wnich  were 
verif ications  of  performaace  reported  in  Ref.  1 or  were 
newly  discovered.  The  uniformity  mentioned  ioove  aiued 
immeasureably  in  the  detection  of  anomalies  since  even  a 
slight  wavefront  variation  could  be  discerned  ia  the  midst 
of  otherwise  uniform  signal  progression. 

Figures  5 and  5 are  two  examples  of  minor 
discontinuities  (affecting  30  to  40  bins  - approximately  100 
kHz)  whicn  appear  to  be  the  result  of  a multi-signal 
amplitude  enhancement  for  a very  short  period  of  time.  This 
effect  was  found  scattered  through  almost  half  of  tae 
thirty-four  data  sets  viewed  via  Anomaly  A.  Figure  7 is  an 
example  of  aperiodic  discontinuities  which  ran  across  almost 
ail  of  the  entice  data  set  bandwidth.  This  change  did  not 
occur  over  the  entire  time  duration  of  the  sat  and,  as 
mentioned  above,  could  not  oe  traced  over  its  entire 
frequency  spectrum.  This  type  of  discontinuity  was  also 
found  in  data  sets  33  and  80.  Both  of  the  discontinuity 
examples  were  anticipated  although  tae  latter  fora  holds  the 
most  promise,  anomaly-wise. 

Fading  was  another  effect  expected  but  tae  examples 
found  were  somewhat  surprising  in  taat  it  was  anticipated 
that  the  entire  freguency  bandwidth  of  a signal,  wavefront 
(890  kHz)  would  oe  affected.  Tnis  proved  not  to  be  the 
case,  with  all  of  the  examples  of  fading  implicating  at  most 
20  freguency  bins  at  a time.  Figure  3 shows  the  signal 
cantered  on  bin  148  fading  in  and  out,  a process  which 
lasted  several  seconds.  • Figure  9 shows  a much  smaller 
time-duration  fading  centered  or.  bin  24. 

A third  type  of  possioLe  anomalous  activity  is  that  of 
tne  shifting,  up  or  down  in  frequency,  of  either  a single 
signal  or  entire  890  kHz  wide  signal  scans.  The  precise 
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Figure  5 - NARROW  DISCONTINUITY  EXAMPLE  #1 
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Figure  7 - WIDE  DISCONTINUITY  EXAMPLE 


FREQUENCY 


Set  33,  Base  Scan  270,  Base  Freq  5 MHz,  Base  Bin  1 


Figure  9 - FADING  EXAMPLE  #2 
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explanation  of  this  phenomenon  is  not  available  oat  it  must 
be  kept  in  mini  that  malfunctions  in  the  receiver  system 
itself  could  be  the  cause  rather  tian  some 
ionosphere-r elated  activity.  Figure  13  shows  a possible 
left  snift  (approximately  3 bins)  and  distortion  of  two 
side-by-side  signals.  Figures  11  and  12  are  examples  of 
apparent  frequency  shifts  of  1 o3  bins  ooth  of  w nich  turned 
out  to  be  data  set  problems.  Since  each  scan  is  composed  of 
two  data  cards  of  binary  coded  amplitudes  (bins  1 - 1o0  an 
card  one  and  bins  161  - 313  on  card  two),  the  order  in  which 
the  two  cards  appear  and,  as  importantly,  the  existence  of 
both  cards  for  each  scan  are  absolutely  vital  to  the  correct 
representation  of  the  signal  wavefront.  For  example,  if 
card  one  of  a scan  is  missing,  tnen  card  one  of  the  next 
scan  is  taken  as  card  two  of  the  former.  Tnis  nas  a domino 
effect  which  will  grapnically  manifest  itself  as  an  apparent 
shift  of  160  bins  in  all  suoseguant  scans.  Figure  11 
illustrates  precisely  tnis  effect.  A manual  search  of  the 
missing  from  the  set.  Figure  12  is  an  interesting  twist  on 
this  problem  ia  that  the  signal  apparently  shifted  160  bins 
to  the  rignt  (or,  down  ia  frequency)  and  then,  taree  scans 
later,  shifted  rignt  bacx  again.  Another  search  of  the 
cards  showed  that  ,as  expected,  two  cards  were  missing,  one 
to  cause  tne  initial  163  Din  rignt  snift  and  a second  to 
return  the  data  cards  to  their  correct  sequence  t.nree  scans 
later.  One  last  example  of  signal  shifting  is  oy  far  the 
most  interesting  of  the  lot.  As  illustrated  oy  Fig.  13, 
there  is  a shift  down  ia  frequency  of  approximately  12  Dins 
(34  kHz)  affecting  tne  entire  318  Dins  (390  kHz  bandwidth) 
of  aper iodicaliy.  spaced  scans.  This  phenomenon  was  seen 
over  the  entire  length  of  the  data  set  at  apparently  random 
intervals.  Data  sets  138  and  143  also  exhibited  tnis 
aperiodic  frequency  shifting,  although  tne  numoec  of  bins  by 
wnach  the  scans  were  sniftad  varied  sligatly  over  the  three 
sets . 
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FREQUENCY 


, Base  Scan  336 , 


Base  Freq  3.2  MHz, 


Base  Bin  100 


Figure  10  - SHIFTING  EXAMPLE  #1 
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Figure  11  - SHIFTING  EXAMPLE  It  2 
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Figure  12  - SHIFTING  EXAMPLE  It  3 
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Figure  13  - SHIFTING  EXAMPLE  //A 


One  final  category  for  possible  signal  anomalies 
consists  of  five  different  phenomena  eaca  of  wain  occured, 
with  one  exception,  only  once  in  the  thirty-foar  data  sets. 
Figure  14  is  an  example  therein  three  different  'anomalies" 
were  present.  (It  shouLd  be  noted  that  ease  scan  = 
"various"  for  Figs.  14  and  15  means  taat  each  of  the  three 
sections  of  the  waveforms  in  the  two  figures  nal  different 
base  scans  and  implies  that  this  in  no  way  detracts  from  the 
usefulness  of  the  figures  as  examples  of  the  affects  in 
guestion.)  Bins  1 througa  40  and  83  through  98  exhibit  a 
very  high  degree  of  signal  attenuation  not  seen  in  normal 
data  sets  (i.e.  data  sets  not  found  to  nave  either  a 
receiver  malfunction  or  test  signal  present).  Bias  41 
tnrough  87  exhibit  normal  signal  characteristics.  Bins  131 
through  200  evidence  a combination  of  unusual  attenuation 
and  normal  signal  activity,  while  bins  231  througi  270  snow 
either  a general  distortion  of  the  signals  present  or  a 
fairly  high  amplitude  level  of  ambient  noise  and/or 
distorted  signals  plus  noise.  None  of  these  phenomena  could 
be  explained  by  receiver  malfunction,  data  card  sequencing 
or  special  signal  characteristics  (i.e.  test  signal,  etc.). 
Figure  15  shows  a data  set  exhibiting  no  continuity  of 
signal  wnatsoevec.  Nona  of  the  signal  "bumps"  can  be 
followed  in  a uniform  manner  as  time  progresses.  rhe  tneory 
that  a serious  receiver  system  malfunction  was  tne  cause  of 
the  distortion  was  borne  out  by  a call  to  tne  cognizant 
authority.  This  malfunction  was  also  verified  for  data  sets 
25  and  54.  The  sawtooth  effect  seen  at  the  middle,  rignt 
hand  side  of  Fig.  16  is  a pattern  not  seen  elsawnere  either 
in  this  set  or  any  of  the  other  thirty-three  sets  analyzed. 
It  shows  wild  fluctuations  in  seven  adjacent  amplitude 
points  (2.8  kHz  apart)  in  one  scan,  whereas  tne  rest  of  the 
scan  appears  to  be  normal.  No  explanation  for  this 
paenomenon  can  be  offered.  Another  unexpected  anomaly  was 
tne  two  flat-topped  signiLs,  centered  at  bins  124  and  12o, 
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igtire  14  - MULTIPLE  ANOMALY  EXAMPLE 


found  in  Fig.  17.  It  appears  to  be  the  result  of  some 
special  attenuation  applied  to  those  frequency  bins  alone. 
Decause,  as  can  be  seen  by  the  rather  "large"  signal 
cantered  on  bin  178,  other  signals  appear  not  to  be 
affected.  The  additional  attenuation  applied  to  this  set 
during  reception  , 18  dB,  is  lower  than  levels  applied  to 
other  sets  which  were  not  similarly  affected.  The  best 
phenomenon  observed  in  tnis  final  category  is  found  in  Fig. 
18.  For  some  reason  (other  than  the  possioility  of  a 
different  form  of  receiver  malfunction  taan  seen  oefore)  an 
otherwise  normal  data  sat  lost  total  signal  continuity 
across  tne  entire  freguency  oandwidth  witn  subsequent  scans 
exhibiting  a hignly  distorted  nature.  while  tnis  could  not 
have  been  caused  oy  a card  sequencing  problem,  its  source 
cannot  be  identified. 


C.  ANOMALY  B ANOMALIES 


The  Anomaly  B prog r a n provided  a different  perspective 
from  whica  to  ooserve  tie  data  sets  since  it  presented  two 
dissimilar  plots  on  tne  screen  at  the  same  time.  The  lower 
traces  consisted  of  five  Anomaly  A style  three-dimensional 
plots,  this  time  only  eigity  points  (or  oins)  wide,  wniie 
tne  upper  five  traces  were  the  two-dimensional  time  vs. 
amplitude  plots  of  Dins  selected  from  the  traces  oelow.  It 
was  felt  tnat  tnis  different  perspective  would  enhance  the 
anomaly  identification  effort,  that  it  did  not  was 
disappointing. 

Of  the  eighteen  data  sets  analyzed  (3,  7,  16,  18,  20, 
21,  33,  3o,  37,  38,  40,  41,  51,  53,  51,  80,  125,  and  139)  , 
Anomaly  3 proved  useful  in  only  two  instances:  (1)  to 
provide  a graphic  presentation  of  the  three  signal  types 
used  in  the  Statistical  Analysis  Section  of  tnis  paper;  and 
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Set  16,  Base  Scan  932,  Base  Freq  9.7  MHz,  Base  Bin  100 


Figure  17  - ATTENUATED  SIGNALS  EXAMPLE 
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Figure  18  - COMPLETE  DISCONTINUITY  EXAMPLE 


(2)  to  illustrate  one  vary  interesting  phenomenon  which 
could  not  be  seen  by  the  Anomaly  A presentation.  As  shown 
in  Fig.  19  this  phenomenon  consists  of  an  undulation  in  the 
amplitude  of  the  five  upper  traces.  This  "lazy  sinusoidal" 
variation  appears  to  affect  some  signals  more  often  than 
others.  The  term  "lazy  sinusoidal"  should  be  use!  with  care 
since  it  is  not  meant  to  imply  taat  there  is  any  strict 
periodicity  to  be  found  in  the  undulation  or  that  the  traces 
look  exactly  like  sinusoids.  However,  many  of  the  signals 
did  evidence  the  more  smooth  (and,  hence,  sinusoidal) 
variation  seen  in  trace  one  in  the  figure. 

Althougn  only  relatively  smootn  waveforms  are  depicted 
in  the  figure,  this  effect  was  identified  over  a large  range 
of  signal  types  including  those  whin  keyed  off  and  on; 
however,  it  was  much  easier  to  see  tnis  phenomenon  on  the 
smooth  signals.  This  variation  was  seen  in  every  data  set 
analyzed  although  it  occured  more  frequently  and  aoticeaoly 
in  some  tnan  in  others. 

The  cause  of  tnis  undulation  could  be;  (1)  variations  in 
transmitter  output  power  levels;  (2)  receiver  system 
processing;  or  (3)  ionosphere-related  phenomena  including 
Faraday  Rotation,  Sporadic-iS  and  fading. 


Set  7,  Base  Scan  588,  Base  Bin  120 


Figure  19  - AMPLITUDE  UNDULATION  EXAMPLE 
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IV. 


STATISTICAL  ANALYSIS 


A.  DISCUSSION 


In  conjunction  with  the  graphical  analysis  approach  to 
identify  signal  anomalies  attnoutable  to  ionospheric 
properties,  an  attempt  was  made  to  glean  additional 
information  iron  the  BRIGHAM  data  by  means  of  statistical 
algorithms;  howaver,  it  should  be  noted  that  no  attempt  was 
made  to  identify  anomalies  by  this  approach.  3 at her,  tna 
programs  MIN/JiAX  1 and  MIN/MAX  2 wara  developed  in  an  effort 
to  determine  now  often  and  by  what  amount  (amplitude  in 
decibels)  both  signal  wavefronts  (scans)  and  individual 
signals  (oins  or  sets  of  bins)  caaage.  It  was  hoped  taat 
this  new  data  mignt  prove  useful  in  the  iavelopamant  of  an 
algorithm  to  adjust  the  myriad  sensitivity  controis  found  on 
the  swept-tuned  receiving  equipment  used  to  acquire  and 
record  the  BRIGHAM  data.  Although  this  effort  nas  nothing 
directly  to  do  with  the  signal  analysis  effort  reported 
previously,  it  did  taka  consiieraule  time  and  effort  to 
write  the  programs  and  craate  tne  dasirad  data  base.  In 
addition,  this  evolution  proved  extremely  haloful  in  tna 
overall  effort  of  understanding  signal  characteristics,  tna 
primary  goal  of  this  study. 

Por  the  purposes  of  tuis  inquiry,  tnree  different  types 
of  signals  (one  that  stayed  "on"  through  most  of  the  signal 
duration,  one  that  stayed  "off"  tnrough  most  of  the  signal 
duration  and  oaa  that  appeared  to  key  on  and  off  throughout 
most  of  the  signal  duration)  ware  cnosen  from  data  sets 


which  were  collected  it  various  times  of  the  day.  The 
quotation  marks  are  used  (here  and  throughout  the  rest  of 
the  paper)  to  imply  a general  trend  noted  over  the  life  of 
the  signal  not  soma  unvarying  signal  characteristic.  Figure 
20  is  an  Anomaly  3 graph  plot  presenting  classic  examples  of 
the  three  signal  types.  Tne  "ON"  signal  (exhibiting  that 
"lazy  sinusoidal"  amplitude  variation  mentioned  in  the 
grapaical  analysis  section)  is  trace  1,  the  "ON/DFP"  signal 
is  trace  2 and  the  "OFF"  signal  is  trace  3.  Traces  4 and  5 
are  merely  copies  of  trace  3 to  fill  out  tne  plot.  A taole 
giving  all  pertinent  data  on  the  twenty-f oar  signals 
analyzed  may  be  found  in  Fig.  21.  In  the  figure,  set 
frequency  refers  to  the  base  frequency  of  the  set  when  the 
recording  was  initiated  while  TDI  means  the  time  of 
interception  for  the  data  set.  An  attempt  was  made  to 
analyze  two  signals  from  each  of  four  different  time 
categories  (interesting  from  the  ionosphere's  standpoint) 

, those  being  sunrise,  nooa,  sunset  and  midnight.  With  one 
exception,  set  53  - 2132LOCAL,  this  effort  was  realized. 


3.  STATISTICAL  ANALYSIS  PaSSRAMS 

The  first  generation  aigorithm  developed,  calLed  MIN/dAX 
1,  provided  for  a scan-by-scan  (down  entire  frequency 
bandwidths  - cross  sections  of  many  signals  - at  fixed 
times)  analysis  of  a variable  number  of  scans  witiin  a given 
set.  dIN/KAX  1 identified  aaca  of  the  peaks  and  valleys 
along  the  scans,  chose  the  hignest  peax  and  lowest  valley 
and  computed  the  average  peak  and  average  valley  heights. 
The  quantization  levels  given  on  the  iata  cards  were  changed 
to  decibels  and  tne  amount  of  attenuation  applied  at  tne 
time  of  acquisition  was  added  back  in  to  give  received 
signal  levels  accurately.  dIN/d  A X 1 proved  to  be  of  only 
small  value  since  signal  wavefronts  vice  signals  themselves 
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TIME 


FREQUENCY 


Set  7,  Base  Scan  200,  Base  Bin  40 


Figure  20  - THREE  MIN/ MAX  SIGNAL  TYPES 
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Figure  21  - MIN/ MAX  SIGNAL  PARAMETERS 
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were  analyzed.  A program  Listing  far  am/MAX  1 may  be  found 
at  tne  end  of  this  paper. 


The  next  generation  effort,  tabbed  MIN/J1AX  2,  focused  on 
specific  signals  (the  three  varieties  delineated  in  the 
discussion  section  above)  and  added  a number  of  more 
valuaDle  statistical  features.  Like  its  predecessor, 
MIN/MAX  2 identified  eaca  of  the  peaks  and  valleys,  chose 
the  highest  peak  and  loMest  valley  and  computed  the  average 
peak  and  average  valley  heights. 

Its  next  feature  Mas  the  determination  of  tMO  different 
types  of  signal  changes:  (1)  the  comparison  of  sample  point 
amplitudes  with  a sliding  rsfarars  point  to  determine 
long-term  changes  in  signal  levels;  and  (2)  the  comparison 
of  adjacent  sample  point  amplitudes  to  determine  short-term 
changes.  The  former  was  initiated  wLth  the  bin  1 poin  u a 5 
tne  reference  point  and  all  suDsegaent  points  compared  with 
it  until  an  amplitude  was  found  which  differed  from  tne 
reference  oy  a previously  agreed  upon  number  of  decibels 
(for  tnis  test,  4,  6,  12  and  13  d3  changes  were  used  in 
succession) . That  point  tnen  became  the  reference  point  and 
the  process  continued  for  the  entire  2.4  minute  length  of 
the  signal.  The  numoer  of  sucn  long-term  changes  and  tne 
elapsed  time  between  each  change  were  recorded.  This  data 
could  then  oe  used  to  attain  the  goal  of  this  effort,  tne 
gauging  of  how  often  and  oy  what  amount  the  various  signal 
types  changed. 

A further  refinement  of  the  procedure  was  the  creation 
of  amplitude  histograms  of  the  elapsed  times  oetween  the 
cnanges  to  illustrate  the  effect  tnat  increasing  the  level 
of  decioel  fluctuation  which  constituted  a change  had  on 
noth  the  frequency  of  cnanges  and  the  distribution  of 
elapsed  times.  Since  the  receiver  sampling  rate  was  25  Hz, 
the  smallest  elapsed  time  possible  between  X dB  fluctuations 
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was  0.04  seconds.  A onef  analysis  of  the  first  set  of 
elapsed  times  led  to  the  division  of  the  time  line 
(abscissa)  into  twenty-one  time  bins  with  the  inclusive 
boundaries  given  in  Fig.  22  and  labeled  SI  through  N21.  As 
a further  refinement,  oins  N1,  N2  and  the  first  segment  of 
N3  were  divided  into  individual  0.04  second  bias  laoeled 
SI-1  througn  N1-7,  N2-1  througn  N2-3  and  S3-1  respectively. 
The  two  histograms  gave  aa  excellent  perspective  from  which 
to  watch  the  changes  in  signal  characteristic  as  the  amount 
of  decibel  fluctuation  tested  was  stepped  from  4 to  6 to  12 
and,  finally,  to  18  d3. 

The  final  output  of  MIS /MAX  2 was,  in  addition  to  the 
peak  and  valley  statistics  mentioned  acove,  a listing  of  tae 
numoer  of  long-term  changes  that  took  place  xn  eaca  of  the 
thirty-seven  time  bins  of  interest,  the  percentage  of  the 
numoer  of  changes  under  each  time  bin  as  compared  with  the 
total  number  of  changes,  (this  allowed  for  a more  valid 
comparison  of  the  over-all  change  uictuce  when  aa  increase 
in  decibel  deviation  reguired  resulted  in  fever  changes 
tabulated  in  the  bins) , the  total  number  of  long-term 
changes  (called  SUM  in  tue  oucput)  and  the  total  number  of 
changes  in  adjacent  sample  points  (tne  short-term  changes 
mentxoned  above).  In  addition,  the  mean  and  standard 
deviation  for  tae  elapsed  times  ioetween  long-term  cnangas 
were  computed.  A listing  for  MIS /SAX  2 may  ce  found  at  the 
end  of  this  paper. 

C.  OBSERVATIONS  AND  CONCLUSIONS 

There  were  a number  of  surprises  which  results!  from  the 
initial  computer  output  £rom  MI N/M AX  2.  First,  the  high 
numoer  of  both  xong-  and  short-term  changes  identified  was 
not  expected.  Out  of  the  possiole  3597  changes  (there  ware 
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HISTOGRAM  TIME  BIN  PARAMETERS 


at  most  that  many  scans,  oc  sample  points,  available  for 
each  signal  of  interest)  tne  number  of  changes  noted  in  some 
cases  reached  as  high  as  2300  for  the  lower  decibel 
fluctuation  levels.  Second,  the  elapsed  time  histograms 
tabulated  showed  a marked  tendency  to  be  heavily  damped 
exponentials  in  nature  witn  the  3.01  second  bin  (N1-1) 
having  the  hignest  number  of  occurences  throughout.  And 
third,  there  was  a high  degree  of  overlap  oetween  the 
statistical  characteristics  generated  foe  the  three  signal 
types  chosen  indicating  that  there  was  no  absolute  way  to 
categorize  a signal  as  "OS",  "ON/OFF"  or  "OFF"  merely  by 
looking  at  its  statistics.  What  was  expected  (and  did 
occur)  was  that,  as  tne  numbec  of  decibels  used  for 
comparison  increased,  tie  exponential  distributions  "flowed" 
to  the  right  indicating  greater  elapsed  times  between  larger 
decibel  fluctuations. 

Figures  23  - 29  represent  tne  bulk  of  the  meaningful 
data  garnered  from  the  SIN/dAX  2 application  to  the 
twenty-four  signals  of  interest.  in  speaking  of  the 
characteristics  which  differentiate  the  three  arbitrarily 
chosen  and  defined  signaL  types.  Figs.  23  - 28  provide  those 
statistics  which  can  be  of  value  in  coming  to  a logical 
conclusion  concerning  tieir  differences.  As  would  be 
expected  from  "ON/OF?"  signals,  tney  snowed  a greater 
difference  between  both  average  peak  and  valley  values  and 
between  "maxmax"  and  "ninnia"  values  than  did  the  other  two, 
implying  more  signal  level  fluctuation.  The  "3F?"  signals 
consistently  snowed  tne  smallest  fluctuations  as  would  be 
expected  from  a series  of  either  highly  decayed  signals  or 
high  but  fluctuating  ambient  noise.  Based  on  these 
statistics  alone  it  woild  be  nearly  impossible  to 
distinguish  with  a hign  degree  of  certainty  between  an  "3N" 
and  an  "ON/OFF"  signal  given  only  the  parameters  in  those 
figures.  Figure  29  provides  both  a summary  of  the  data 
concerning  tne  number  of  changes  tabulated  for  each  decibel 
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Figure  27  - "ON/OFF"  SIGNAL  ELAPSED  TIME  STATISTICS 
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Figure  28  - "OFF"  SIGNAL  ELAPSED  TIME  STATISTICS 
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comparison  level  and  a general  means  of  distinguishing 
between  the  three  types  of  signals.  Althoagi  there  are 
definite  differences  between  the  average  figures  for  the 
various  dS  levels  and  some  differences  in  the  ranges  of 
values,  the  overlap  encountered  in  the  ranges  suss  signal 
differentiation  by  this  means  risky.  In  snort,  wnat 
distinguished  the  "OH"  from  me  "ON/ OP?"  and  "OP?"  signals 
was  thoroughly  subjective  and  only  loosely  corroborated  dv 
the  statistics  generated. 

Additional  information  can  be  retreivea  from  Figs.  2b  - 
29  wnich  provide  the  numoer  of  cnaages  (both  Long-  and 
short-term)  found  waen  tne  decioel  differential  was  "walked" 
from  4 to  13  dfl  and  the  percentage  of  long-term  changes 
which  fell  into  the  most  meaningful  of  tne  tnrrty-seven 
histogram  bins  (for  elapsed  times  between  long-term 
changes)  . As  is  clearly  shown,  the  number  cf  cianges,  n'JM 
and  ADJ,  decreased  as  the  dB  level  increased  while  the 
percentage  figures  in  coiuans  N1  and  N1-1  evidenced  a 
movement  away  from  the  first  elapsed  time  bin.  Both  of 
these  character istics  imply  that  tne  average  time  between 
signal  level  changes  at  tne  various  d3  settings  increases  as 
tne  number  of  d3  increase.  To  be  specific,  Dy  taking  tne 
data  from  Fig.  30  it  couLd  ae  calculated  that  tne  average 
time  elapsed  oetween  4 and  6 d5  fluctuations  foe  the  tnrea 
signal  types  were: 


Signal  Type 

4 dB 

3 

dB 

ON 

. 195  sec 

. 3 o 2 

sec 

ON/OFF 

.033  sec 

.132 

sec 

OFF 

.110  sec 

. 229 

sec 

5o 


— 


SET 

8 I N 

NCB 

MEAN 

ALPHA 

CHISC 

DF 

VARIANCE 

3 

136 

4 

0.116 

_ 

__ 

2 

136 

6 

0.226 

— 

— 

— 

— 

5 

154 

4 

0.067 

— 

— 

— 

— 

3 

154 

6 

0.108 

0.93 

1.43 

1 

C.  012 

3 

169 

4 

0.088 

0.87 

0 .004 

1 

0.009 

3 

169 

6 

0.  179 

0.34 

9.18 

5 

G.  C94 

7 

8 

4 

0.102 

— 

- 

— 

— 

7 

8 

6 

0.279 

0.18 

5.88 

5 

C.432 

7 

81 

4 

0.072 

— 

— 

— 

— 

7 

81 

6 

0.106 

— 

— 

— 

— 

7 

91 

4 

0.229 

C .34 

3.26 

5 

0.154 

7 

91 

6 

0.985 

0.18 

2.73 

5 

5.39C 

20 

26  5 

4 

0.096 

— 

— 

— 

_ 

20 

265 

6 

0.171 

0.16 

8.35 

5 

C.  183 

2C 

238 

4 

0.067 

— 

— 

- 

— 

20 

238 

6 

0.095 

— 

— 

— 

— 

20 

259 

4 

0.090 

0.94 

0.34 

1 

0.009 

20 

259 

6 

0.189 

0.50 

4 .CO 

5 

0.07  1 

21 

150 

4 

0.  100 

— 

— 

— 

— 

2 1 

150 

6 

0.157 

C.C45 

7.36 

5 

C.  548 

21 

273 

4 

0.091 

0.85 

1.52 

1 

C.01C 

21 

273 

6 

0.153 

0.49 

4.21 

5 

0.048 

21 

187 

4 

0.067 

— 

— 

— 

— 

21 

167 

6 

0.143 

0.44 

2.33 

5 

C.046 

53 

19 

4 

0.436 

— 

— 

— 

— 

53 

19 

6 

0.748 

— 

— 

— 

— 

53 

8 

4 

0.CS5 

— 

— 

— 

— 

53 

a 

6 

0.136 

0.065 

12. 1C 

C.284 

53 

50 

4 

0.085 

1.05 

0.18 

l 

0.007 

53 

50 

6 

0.174 

C.28 

9.32 

5 

C.  108 

EC 

82 

4 

0.162 

0.18 

10.65 

5 

C.  146 

EG 

82 

6 

0.411 

— 

— 

— 

— 

80 

60 

4 

0.101 

— 

— 

- 

— 

EC 

60 

6 

0.218 

0.18 

14.58 

5 

0.264 

60 

73 

4 

0.081 

1.05 

0 . 6 C 

1 

0.006 

60 

73 

6 

0.166 

— 

— 

— 

— 

£1 

4 

4 

0.074 

— 

- 

• 

— 

61 

4 

6 

0.106 

— 

— 

— 

— 

6 1 

72 

4 

0.062 

— 

— 

- 

— 

81 

72 

6 

0.087 

— 

— 

— 

— 

61 

11 

4 

0.C69 

— 

— 

— 

— 

61 

1 1 

6 

0. 113 

0.76 

2.70 

5 

0.017 

126 

57 

4 

0.424 

0.72 

1 .16 

5 

C.  250 

126 

57 

6 

0.79  5 

0.62 

6 .50 

5 

1.019 

1 26 

74 

4 

0.111 

— 

— 

— 

— 

126 

74 

6 

0.245 

0.28 

10.38 

5 

0.214 

126 

8 

4 

0.152 

— 

— 

— 

— 

126 

8 

6 

0.452 

0.51 

2.68 

5 

C.40C 

KEY  : 

NCB  = 

MI N/M AX 

AMPL ITUCE 

CHANGE 

VAR  IAELE 

ALPHA  = VALUE  CF  AlPHA  AT  MIN  CHI  SQUARED 
CHISQ  = MINIMUM  CHI-SQUARED  VALUE 
OF  = DEGREE  CF  FREEDOM  FC«  CHISQ 


Figure  3 0 


CURVE-FITTING  DATA 


This  data  confirms  taat  tha  "ON/OFF"  signals  fluctuated 
aore  rapidly  than  tha  otner  two  categories  wnile  the  "ON" 
signals  proved  to  ba  more  stable  than  the  "OFF" . 

These  results  imply  that,  if  tae  agency  that  operates 
the  swept-tuned  receiving  system  used  to  acquire  and  record 
tne  BRIGHAM  data  desires  to  acquire  FSK-tyoe  signals 
(roughly,  the  "ON"  variety),  it  could  get  away  with  changing 
tne  sensitivity  adjustments  much  less  frequently  than  if 
Morse-type  signals  (rougaiy,  the  "ON/OFF"  variety)  were 
desired.  This,  plus  tne  rough  data  provided  in  Fig.  23, 
could  provide  a good  data  base  from  which  a dial  controlling 
algorithm  could  be  adjusted  to  meet  the  needs  of  the  signal 
types  of  interest  at  any  time. 

As  was  mentioned  above,  MIN/MAX  2 provided  for  the 
calculation  of  the  mean  and  standard  deviation  of  the 
elapsed  times  tabulated  for  eaca  of  forty-eignt  signals  (the 
twenty-four  signals  of  interest  at  noth  '*  and  6 d3 
comparison  levels).  Altiough  tha  mean  values  were  accurate 
the  standard  deviation  values  were  not  since  tae  formula 
used  in  their  calculation  is  valid  for  only  normal  (or 
Gaussian)  distributions.  Since  the  elapsed  time  histograms 
were  apparently  exponential  in  nature,  anotner  form  of 
variance  determination  was  required.  This  effort  called  for 
a curve-fitting  technique  to  identify  tne  histograms  as 
being  of  a particular  distribution  type,  use  of  tae 
chi-squared  test  for  goodi ess-of -f it  and  manipulation  of  the 
cnosen  distribution’s  meaa  and  variance  formulas  to  derive 
the  information  desired.  Attempts  were  made  to  fit  three 
distribution  functions  to  the  histograms  (Gamma,  Poisson  and 
Geometric)  with  the  Gamma  function  proving  to  be  the  most 
successful.  The  Gamma  uistrioution  function,  with 
parameters  3</3  < is  defined  as  follows: 
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f ^ " 1 £*-1  "x//9 

([»/T>  x a p,  x>D 

f(x)  = < 

I 0,  xO 

Where,  Mean  = c*/3  , Variance  =fr/2>  . 

An  algorithm  was  written  wnich  real  in  the  observed 
histogram  frequencies,  calculated  the  expected  Gamma 

fregaencies  (using  an  Intarnational  Mathematical  Subroutine 
Library  subroutine  caLled  MDG&M)  and  calculated  the 
chi- squared  value  (essentially,  the  sum  of  the  least  squared 
differences  between  taa  ooserved  and  expected  frequencies) 
for  use  in  the  goodness- of- fit  test.  Ihe  chi-sguared  value 
was  then  compared  wita  the  appropriate  value  found  in  tne 
chi-sguared  distribution  taole  (appropriate  as  to  degrees  of 
freedom  and  significance  level);  if  it  fell  below  that  table 
entry,  the  fit  was  considered  good,  if  not  the  fit  was  bad. 
As  discussed  in  Refs.  7 and  6,  a good  rule  of  thumb  is  that 
the  smaller  the  chi-sguared  value  for  a given  number  of 
degrees  of  freedom,  the  better  tne  fit.  As  a word  of 
caution,  even  an  excellent  fit  (extremely  small  chi-squared 
value)  does  not  guarantee  tnat  the  distribution  is  really 
the  one  being  tested,  it  necely  states  tnat  the  data  do  not 
present  sufficient  evidence  to  contradict  tne  hypotnesis 
that  the  nistogtams  possess  Gamma  distributions. 

For  the  curve- fitting  at  hand,  knowledge  of  tie  mean  and 
an  iterative  variation  in  c*  permitted  an  iterative 
calculation  of  /!>  thus  providing  all  of  the  necessary 
information  for  tie  algoritna  discussed  above.  Tie  process 
was  refined  until  tne  saalLest  chi-sguared  value  could  be 
identified.  It  saould  oe  noted  that  only  the  Ml  through  M15 
time  bins  were  used  since  very  few  elapsed  times  ever  fell 
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beyond  five  seconds  and  that,  as  the  individual  histogram 
composition  dictated,  either  seven  or  three  bins  were 
actually  used  in  the  algorithm  (following  a suggestion  that 
bins  containing  five  or  fewer  counts  be  combine!  to  enhance 
the  accuracy  of  the  test) . The  bin  limits  implied  that  the 
degrees  of  freedom  used  were  five  or  one  respectively. 

Of  the  forty-eight  Histograms  to  be  fitte!,  six  had 
counts  in  too  few  bins  to  be  useful  (would  have  inpliea  zero 
degrees  of  freedom  - not  found  in  tne  taole) , eigiteen  were 
nad  fits  and  the  remaining  twenty-four  provided  good  fits  to 
the  Gamma  distribution.  Figure  30  gives  the  results  of  the 
curve-fitting  effort.  Knowledge  of  rx.  and  tne  mean  implies 
Knowledge  of  (3  which,  ia  turn,  implies  (by  definition) 
knowledge  of  the  variance.  A combination  of  tie  mean  and 
variance  values  saould  be  useful  in  the  determination  of  how 
often  the  sensitivity  dials  on  the  receiving  system  should 
be  adjusted  for  tae  desired  signals  to  be  acquired. 

Figure  31  provides  a breakdown  of  the  <X  values  vs. 
signal  type  and  amplitude  change  variable  for  the 
twenty-four  histograms  which  "fit"  the  Gamma  distribution. 
It  was  felt  that  tais  would  permit  the  grouping  of 
Histograms  with  respect  to  some  combination  of  the  three 
categories  amove.  The  six  ranges  snown  in  the  figure  impLy 
that  any  grouping  by  alone  would  be  tenuous  indeed  as 
they  overlap  quite  heavily.  The  only  discernaole  trend  is 
tne  apparent  lowering  of  the  avecage  <x  figure  as  the 
amplitude  change  variable  increases  from  4 to  S IB.  This 
would  seem  to  contradict  the  defining  equations  for  the 
Gamma  distribution  whicn  imply  tnat  a lower  cx.  implies  a 
more  "daaped-ex ponen tial ’ fora  for  the  distribution  (when, 
in  fact,  the  more  damped  case  is  found  at  the  4 d3  change 
variaoie) . however,  the  shape  of  the  distribution  is  also  a 
function  of  (3  wnich,  as  iefir.ed  in  'lean  = , tends  to 
affect  the  distribution  snape  in  a more  radical  fashion 
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than  does  C<  . in  fact,  as  tae  change  variable  increases, 
the  mean  increases  (see  Fig.  30)  and  c<  seems  to  decrease 
(see  Fig.  31)  implying  a vary  pronounced  increase  in 
which,  in  the  defining  equations,  is  found  to  be  most 
prominent  in  the  expoaeatial  term;  jut,  this  implies  a 
distribution  shape  lass  heavily  dampled  in  the  lower  (4  dB) 
change  variable  level.  Tais  trend  of  decreasing  values, 
then,  was  to  be  expected  and  appears  to  be  the  only 
meaningful  product  of  the  analysis  of  the  Fig.  31  data. 

is  an  aside  an  the  suojact  of  signal  continuity  during 
the  brief  recorded  time  span  of  toe  waveforms  (around  2.4 
minutes) , it  should  be  noted  that  data  set  nuaoer  7 bins  91 
("ON/OFF")  and  91  ('OFF")  exnibited  vary  unusual 

characteristics  when  compared  with  the  other  signals  in 
tneir  classes.  As  snown  on  Fig.  27,  oin  81  stiLl  had  over 
1000  changes  at  13  dB  difference  while  tae  otner  "OH /OFF" 
signals  had  fallen  off  to  much  smaller  numbers.  In  fact, 
that  oin  had  two  changes  at  54  dB  wnen  ever-increasing 
numbers  were  applied.  The  "OFF"  signal,  bin  91,  evidenced  a 
similar  degree  of  consistency  as  it  nad  many  fewer  4 dB 
changes  than  the  other  "OFF"  signals  wnile  falliig  off  much 
more  guickly  as  tne  number  of  decibels  increased  (see  Fig. 
28) . In  fact,  this  signal  nad  only  a few  more  than  1700  2 
dB  changes  implying  exceptional  amplitude  continuity  when 
compared  with  any  of  the  twenty-three  other  signals 
analyzed.  To  give  a more  grapaic  representation  of  tnis 
phenomenon,  Figs.  32  and  33  compare  the  long-term  signal 
changes  found  in  data  sets  3 and  7 with  the  latter  proving 
to  De  the  mot  e consistent  of  tne  two,  although  even  set  3 
nas  a good  measure  of  uniformity  througnout  its  3600  scans. 
The  base  scans  for  each  of  the  five  scan  sets  (whica  show 
the  beginning,  middle  and  end  of  each  data  set  - tirae-wisei 
provided  at  tne  rig at  of  the  pictures. 


o 2 


are 


FREQUENCY 


Base  Freq  11  MHz,  Base  Bin  1 


Figure  33  - LONG-TERM  SIGNAL  CONTINUITY  - SET  7 
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V.  CDUIlESrS  AND  REC3HHEN  DATIDNS 


3f  the  two  signal  anomaly  identification  techniques. 
Anomaly  A proved  to  be  the  most  fruitful.  Narrow  and  wide 
wavefront  discontinuities,  fading,  frequency  shifting  of 
entire  wavefronts,  frequency  shifting  of  individual  signals 
and  several  miscellaneous  varieties  ware  seen.  Two  of  the 
above  (discontinuities  and  shifting)  were  reported  on  in 
Ref.  1 with  the  rest  being  unique  to  this  study.  several 
examples  of  pseudo- anomalies  were  included  to  perform  tne 
dual  function  of  warning  against  tae  false  assumption  of 
anomalous  performance  simply  because  the  plot  seemed  to 
imply  it  and  to  provide  pictures  witn  tne  narrative  in  such 
cases.  Analysis  with  Anomaly  B proved  to  be  relatively 
unproductive  as  only  the  interesting  and  frequently  seen 
pnenomanon  of  signal  amplitude  undulation  (reported  in  Ref. 
1 and  expanded  upon  nerai  - oeing  noted.  The  diverse  and 
amazingly  complex  nature  of  tne  ionospnere  and  its 
interaction  with  its  environment  precludes  tie  precise 
assignment  of  a cause  to  every  effect  (i.e.,  anomaly). 
Suffice  it  to  say  that,  etcept  for  those  anomalies  wnicn  can 
be  explained  away  cy  equipment  malfunction,  notning  has  been 
found  in  this  effort  wnicn  will  simplify  tne  modeL  of  the 
ionosphere  as  an  ever- changing  medium  in  wnicn  literally 
anything  can  happen  to  a radio  signal  wnich  traverses  it. 

The  statistical  anaLysis  worx  performed,  from  creation 
of  tne  MIN/MAX  algorithms  to  the  tedious  c ur ve- f it t i p.g 
effort,  resulted  in  a wealth  of  predominantly  raw  data 
wnicn,  witn  the  mean  and  variance  values,  should  provide  a 
foundation  upon  wnicn  a receiver  system  sensitivity 
dial-controlling  algoritam  could  be  refined. 

o5 


E t was  noted 


in  the  curve-fitting  section  that  it  is  perhaps  tae  weakest 
link  in  the  statistical  study  since  the  chi-sguared  test 
guarantees  neither  that  the  elapsed  time  a ist o gr a ms  ate 
Gamma  distributions  noc  that  there  Bight  exist  other 
distributions  which  wouLi  "fit"  the  histogram  data  as  well 
as  the  Gamma  function  did.  How  e vac , a statistically 
reasonable  argument  can  be  made  that  the  histograms  (at 
least  half  of  them)  are  probaoly  Gamma  in  nature  and,  hence, 
the  variance  values  calculated  can  be  accepted  as  accurate. 

Should  additional  wort  oe  done  with  the  BRIGHAM  data, 
the  following  remarks/ recommendations  could  prove  useful: 

(1)  all  of  the  compLete  data  sets  in  hand  nave  been 
thoroughly  anaLyzed  wit  a program  Anomaly  A - additional 
effort  along  thus  line  uijit  oe  a waste  of  time  and  energy; 

(2)  althouga  nearly  twenty  of  tae  data  sets  in  hand  have 
not  been  analyzed  with  program  Anomaly  3,  the  disappointing 
results  thusf ar  mignt  act  as  a warning  as  to  expected  future 
rewards  from  additional  effort; 


(3)  if  Anomaly  A and/or  Anomaly  3 ace  to  oe  jsad  again, 
the  bAKSCN  feature  found  in  the  dErASYdoDL  subroutine  (it 
provides  for  tne  backing  over  of  a specifiai  number  of 
files,  or  data  sets)  should  be  nada  to  worx  as  a time/lacor 
saving  device  - the  proolem  is  in  tne  magnetic  tape  unit, 
not  the  program; 


(4)  the  recommendation  (made  initially  in  Ref.  1)  for 
the  creation  of  a computer  algoritnm  to  seaccn  out  the 
anomalies  in  an  automatic  and,  hence,  more  rapid  manner  tnan 
the  visual  approach  is  guardedly  endorsed  herein  - although 
the  speed  at  wuich  the  data  sets  could  oe  analyzed  would  be 
improved  by  at  Least  a.n  order  of  magnitude,  tne  extreme 
complexity  of  trying  to  tell  the  computer  what  to  look  for 


as  well  as  the  loss  of  tnat  very  intangible  but  valuanle 
human  discretionary  power  could  combine  to  made  this  effort 
a quagmire; 

(5)  should  additional  work  oe  done  utilizing  the  HIN/HAX 
routines.  Anomaly  3 shouli  oe  used  first  to  enhance  the 
accuracy  in  selecting  tie  different  signal  types  to  oe 
analyzed  - that  this  was  not  done  ia  the  present  study  aay 
have  resulted  in  the  rather  large  ranges  found  in  the  date; 

(6)  the  accuracy  of  any  future  c ur ve- f itting  work  couid 
be  enhanced  by  the  inclusion  of  the  tfeibull  and  Lognormal 
distributions  (or  aay  of  the  1 e v el- cr ossi ng  distributions) 
to  the  list  of  possiole  curves  to  which  the  eLapsed  time 
histograms  are  matched  because  sets  of  such  data  can  be 
fitted  to  more  than  one  iistrioution  with  equal  accuracy  and 
validity  - this  was  not  done  here  oecause  the  suggested 
distrioutions  were  compLex  eaouga  and  the  time  ren  aiming 
short  enough  that  a good  effort  could  not  have  oeen  put 
forth. 


1.  S^oradic^E:  thin  horizontal  layer  (or  patcm  of  hign 

electron  density  subsided  in  the  regular  E layer 
resulting  in  signal  deflection  and/or  absocption. 

2.  Night^bay  Asymmetry:  Less  attenuation  at  aigit  than  day 

due  to  difference  in  the  ionizing  effect  of  the  sun 
resulting  in  discontinuities  during  the  suar ise/sunset 
periods . 

3.  Seoaa  cine  tic  Non-Re cigcocity : different  attenuation 

effects  on  signals  travelling  west  to  east  rather  than 
east  to  west  due  to  tne  earth's  magnetic  field  lines. 

4.  E.2ii!L  Cat  Absorption  L1Zh.L:  aosorption  of  radio  signals 

passing  tne  polar  regions  - due  to  a conoir.ation  of 
solar  flare  activity  and  the  shape  of  tne  earth's 
magnetic  field  lines. 

5.  Auroral  2i§2iiy:  origatiy  colored  aortaera  latitude 

phenomena  disruptive  to  communications  - also  as  the 
result  of  a combination  of  solar  flare  activity  and  the 
earth's  magnetic  field  lines. 

o.  Faraday  Rotation:  any  Linearly  polarized  wave  traveiiina 

in  the  direction  of  a magnetic  filed  results  in  its  two 
circularly  polarized  components  travelling  at  different 
velocities  and  tnus  the  plans  of  polarization  will 
rotate  witn  distance. 
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if  a radio  w a/a  comas  into  contact  with  a region 
of  inhomogeneous  rafractiva  indices  the  wave  could  be 
trapped  between  two  layers  and  guided,  as  in  a leaky 
wave  guide,  away  fron  its  intended  dastination. 

fLisliaa  IDallinjer  FfractL:  sadden  ionospheric 

disturbance  (S.  I.  ) which  produces  a complete  radio 
"fade  out"  lasting  from  a few  minutes  to  an  hour  or 
more  - caused  Dy  solar  flare  activity. 
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AM.MALY  A - PROGRAM  1C  DRAW  A TP  R t £ DlKcNS  ICRAL  PICT  (TIME, 

FREcL  £Nv,Y , AMPLITUDE)  GF  HF  SIGNALS  PAVING  TRAVERSED  TPt  ICNCSPPERE 
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PLOT  ROUTINE 


aa 


. c 


*de  f i ne 

XINT 

320 

#de  f $ n« 

YINT 

6000 

*def i ne 

OELY 

3200 

#de  f i ne 

OEL  X 

16 

*de  f i ne 

AMPY 

6000 

*de  f i ne 

NLINE 

1 

#de  f i ne 

PI 

3. 

•define 

NP1 

ISO 

int  zero  0; 

int  one 

t ; 

st  rue  t ( 

int 

not  ! 

int 

line; 

i nt 

Ids; 

i nt 

1 dy  > 

i nt 

cot  r ; 

int 

nsot  ; 

i nt 

soot l ; 

int 

soot  2 ; 

i nt 

spot  3; 

i nt 

h r > m i n . 

sec  f 

i nt 

nf  t ; 

i nt 

1 aa  l 

int 

srl 

i nt 

i d I 3)  J 

i nt 

i d; 

i nt 

sc  1 f 

int 

no  r day  t 

yr} 

> head. 

»h ; 

int  hbut(300); 
s t rue  t data  ( 
I nt  » ; 

Int  yd! 

) ; 

st  rue t dat  a *d» 
i n t dbo  f 1 «00  01  ! 

I n t np  J 
int  nip  9 ; 


Int  i dew  , pdev  .se1*»selv»Mas»»tdev; 

»« i n ( arac  t bpov ) 
int  »*aPQV« 

I 

int  i » j » n ; 
ebar  *c  s : 

t Mildew  — ooan  t * /daw  / spd"  . I ) ) < 0)( 
PC i nt f ( "e snnot  open  spp  »"); 
esitO; 

> 

iMfndev  = open  ( " /di»v  / r vp  ” , 1 1 ) < 0 ) ( 
or i nt f ( "c annpt  open  rvn  «"1>" 
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esitm 

1 

i f ( ( i de w * open ( " /dev/rmt 6" , 0 ) ) <0 ) < 
print  t ( "cannot  open  rmtb  »")? 
e.itO? 

» 

i f (eroc  > t ) I 

c s = arov  (11; 

n : 0 J 

while  ( «es  > - 'O'  tS  *C*  « = '9') 
n 2 n • 10  ♦ *C$M  - 'O'; 
n 2 » III 

for  ( i sOl  i <n ? i ♦ ♦ ) 

inn(idev,  dbu(«O00)? 

print  f (“number  o*  records  skipped,  Xd»",n); 

I 

while((tnp{idev,hbuf»000))»0)( 
h s Hbuf? 

1 = o; 

lord  s<);  I < h ->  line!  (♦♦)( 
r.  s t no ( i dev , ftdbu H j I • 000)? 

I st  ft ( h- >no t ) ? 

» 

sc1«  = 010;  setv  = 05? 
bi as*=50; 

no  2 h->not  * h->line? 
sc  a I e ( ) ? 
p I ot ( ) ? 

C ver s ( odev , 020 ) ? 
s 1 1 v ( odev , tone ) ? 

I 


*< l€  * \ ne 

NHYT  2b« 

*de  f < ne 

NHL  0 0 

•define 

NSL  1250 

*de  f < 

D9AW  1 

St  rucl  lot ( 

i nt 

r • 

C h«p 

**pb? 

vine? 

c *ar 

«di  r? 

1 nf 

yfl 

< nt 

(ml 

♦ nt 

* * 1 Ink? 

t nt 

*b 1 ink? 

I i t ab  i 7*MRV  r ) , • in,  • i a? 
ebar  nb  INB t T 1 ? 


1 


4 nt  *dp (2000) ; 


P I ot ( ) 

( 


struct  (pt  »s« 
i nt  i . j r i e ; 
sort ( ) ; 

for(i  = 0;  i < NBYT;  iff) 

pb ( < 1 = o; 

f or ( 4 = 0;  4 < NHIK;  iff) 

wr4tefpdev»nbrNBYT); 
ip  - OS  <a  = 4 1 ab ; 
i = OS 

ford  = NSIS  i > -IS  i — )< 
whilefj  < no  Mt  • rip  f j 1 
< e - s 4 P ( dp  ( j ) ) : 

4 < ( 4 A ==  n 


//clear  plot  buffer 


//move  to  too  cf  plot  area 
//pint  scan  line 


= = 4)4 

//set  up  plot  point  for  4 n t e rpo 1 a t 4 on 


ret  urn; 

> 

n)b(i);  //set  uo  plot  buffer 

wr 4 t e (ndev r ob » NR Y T ) ; / /p I o t line 
t 
V 


4 nt  | sort  # 

sort ( ) 

4 

reoi s t er  4 , k , t S 

4 -OS 

forirf  = dbu  f S d < tdbuffnp*2l;  dff) 
dp  lit*)  = Ud  ->  vd); 

k = np ; 

while  4 k = >>  1 ) 4 
i sort ; 

while  4 j sort  ) 4 
Jsort  1 0? 

for  4 4=0;  i«  (nprk);  4 ♦ ♦ ) 

4 f ( *do  I i ) s •dolitkl  )4 
t=dolil  ; 
dp  4 1 ) =do  1 4 fk)  ; 
dp  4 4 f k I = t « 
j sort ♦ ♦ ; 
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s i p (dDi  ) 

i nt  *dp i i 

( 

in t i t i e • 
s = dp  i I 

\ c * 2 1 ? 

yi  2 *dpi**»  *i  2 *dni# 

if(*4a  < R (dbu f (2*no)  ) ) < 

* r - * s ♦ ♦ ; 
y r 2 ' * s S 
i f (y r & DR Art)  < 

1 f ( ( i nc  2 v i “ y r ) > = 0)  ( 
i f ( ( i nc y 2 i nc ) 22  0 ) 
di r w 2 nip; 
e 1 se  { 

f or ( i 2 0;  ((incy  2 (inc/(nio  >>  i))) 
dir*  2 ( 1 <<  i)# 

) 

f y = y r ? 
c*  2 0?00? 
i c * 2 o ; 

ie  2 st^cMvl-incy»incy»fy#xi»d<rx»cx)f 
i f ( i 22  1 ) 

rpt  urn ( 2 ) ! 


i f ( --do i >2  dhu  f ) ( 
i f (y i R DM  Ad ) ( 

y \ 2 *do i •• I 
x 1 2 *dP i ! 

i f ( ( i nc  2 v*  • v 1 ) >=  0)1 
i f ( ( i nc  v 2 i nc ) 22  0 ) 

dir*  2 -n i p ! 
d \ se  < 

for(i  =0;  ((incy  2 ( i nc / (n i d >>  i)))  22  0);  t ♦ ♦ ) i 
dir*  2 -(I  <<  i ) * 

) 

f v 2 y I i 
c * 2 1c** 

id  2 staefc(y<*incv*incy#fy#*i-l#dirxrcx); 

i f ( i € 22  1 ) 

rdf  urn ( I ) ; 


rdf  urn ( 0 ) 1 


stacir(a»b»c*dd»d»  f) 

inf  a#h»c»dd  #<*»♦; 


< 

Int  »s; 
int  i ; 

struct  iot  ♦ *; 

la  ->  y - a ; 
i a * > vine  = b; 
i a - > y f r c i 
i f (dd  >=  MBIT  I IS  dd  < 0 ) { 

printf(”bad  scale  « s*di*  " , dd ) I 
e * i t ( ) ; 


) 

I a 
I a 


»ob  - dd  ■ 
->  xdl r = e; 


ob ; 


ia  •»  exp 


t; 


I f (ip  OH 
io=i t ab; 

ip->fllnk=in->blinksO; 

fata; 

i a->b 1 (nk=io; 
i a - > f I InksO; 

) 

el  se< 

s=la-»bl Ink; 

S*>fl ink:ij! 
t f( ia-»f| ink  ::  0)1 
s - la; 

I f ( ♦ ♦ i a > = Hi  tab (7*NBYT)  ) ( 

prlntff’itab  overflow  a " ) ; 
return(l); 

> 

i a->bl I n k = s ; 
i a->  * 1 i nV  = 0; 

» 

e I se  ( 

s = i a ->  f link ; 
s *>  blink  s la! 
ia  ->  flink  s 0; 

I a = s ; 

» 

) 


n i b ( s I ) 
I nt 

( 

I nt 
i nt 
whi 


a I ; 

I t i I n; 

* s ? s = in! 

I e ( s ) ( 

i f ( s ->  y f <0 ) 

I f ( ( s=f ree ( s ) I = s 
return; 

*s  ->  *pb  =s  ->  c «r; 


0) 


I f (si  3 • 

i = ( (n 


v’  I 

3 - > 


xd i r ) > 0 ? 


//setup  nlottinq  buffi 
//delete  point 


//line  break 
i : -n ) ; 
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1 


//left 


i f < n < 0 ) 

fori  i = o;  i < i ; J to  { 
i f ( s ->  c «o  " OH 

*9  ->  mb  = * If 

9 •»  c xo  - If 

» 

»s  ->  xob  -i  (s  ->  cxn  :<<  I)f 


else 


for ( ) = o;  i < if  »♦»)( 

>s->mh  s!  (s->c«p  ->>  I)! 

» 

s->y  = - s->v i nc  f 


> 

) f ( s I « = s->y  f ) 

9->v  f - "If 

9 - S">  f H nk  f 


/ /end  of  PO  i n t 


f reel s ) 

lot  *s  f 

< 

int  * t f 
lot  If 

struct  lot  *zS 

»s  •»  xpb  s Of 
I f (9-»bl  ink  OH 

ip  = s">  f I i nk  f 
i p->b 1 i nk  — Of 
t = i Df 

> 

e I se  < 

t = 9->blinkf 
t">flink  - 9*>f I inkf 

t = 9 -»  flinkf 
i f ( t ==  OH 

9 •>  flink  = i a f 
i a - s f 
ret  urn ( t ) f 
I 

t •>  blink  - 9 ->  hlinkf 

I 

S • > f I ink  - iaf 
!*>bl Ink  2 la->blinkf 
i a = 9 f 
ret  urnl t ) f 
I 


sc  a I e ( ) 

< 

st  rue  t data  »s  f 
int 

int  «i  * I £ 0 I » dy ( ^ A I * 
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s 2 dbu f 5 min* 

r m i ny  = 

077777; 

d=s ; 

for(i  = 0;  i < 

np  S i ♦ ♦ ) ( 

min*  = ( (c 

= s - > * ) 

< min*?  c ; min*); 

m i ny  = ( (c 

= S •>  yd)  < miny?  c : miny) 

s*  ♦ ; 

) 

i f ( m i n * > 0 ) 
min*  ~ 0! 
if(miny  > 0) 
m i ny  = 0 ; 

S=dbu  f ? 

f op ( 1 = 0 ? i < 

np ; i ♦ f ) ( 

dm  r s -> 

yd  K 1 ? 

s ->  * s- 

min*# 

s •>  yd  = • 

miny; 

s ->  * =/ 

SC  1 * i 

s ->*  = ♦ bias*# 

it(5->.  >= 

NHVT ) ( 

print  M"ov*rf  low 

nb  »"  i; 

e«1 t ( ) 

1 

s ->  yd  -/ 

0 

sc  1 y ; 

1 * ( s->yd  > 

NSl) 

s->yd 

• NSl-1  ; 

s -»  yr)  = t 

017777b; 

s -»  yd  =! 

dm ; 

s ♦ ♦ ; 
i 

> 

( np  ( idf,btjf>nbvt») 

Int  1 d f § *bu f t nbv r r ; 

< 

i nt  1 * t » n, c i 
Struct ( 

char  el«c?.c?.cu; 

> c < 1 12001 , *s; 
s = c f ; 

if((n  s Hf ,c  t »rbvt« U > 0){ 

torfl  s 0 l i < nb  V t (•/  <1 ; i**M 
c = s -»  tZ  <<  ?; 
t = c < < 10; 
t =!  s ->  c 5 <<  6; 
t -\  s ->  ct; 

*♦  ♦ ; 

bu  1 1 f l = t ; 

} 

) 

> 

char  »cbp*  * 1 oe » t M I 521  » ch  1 1 0)  ; 


c onv ( v a 1 ) 

i nt  vat; 
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< 


I n t a 7 

i Ha  = val / 1 0) 
conv (a ) < 

*cbo*  ♦ = val  t 10  ) 'O'; 
I 


cone (c  1 » c <? » n ) 
chap  »c?) 
i n t n, >c  I ; 

< 

i nt  I in; 

lord  = o;  i < n;  I 4 4 ) ( 

I Mi  i=  0) 

" = ( * c I l 07700)  >>  h; 

else 

">  = *c  I 3 007  7; 
ill™  :=  01?) 

*c<?  - 'O': 

lf(m  >=  0 1 33  7i  <:  Oil) 

*c2  = ' 1 * ♦ 7.  - | ; 

i Me  >=  0?1  &3  7i  < z Oil) 

»c2  s ’ A • 4 rn  - 0? I ; 
i Mm  >=  00  | JR  71  <:  OS  l ) 

* c ? = 'J'  4 7i  - 0 0 1 ; 

i f (ti  >s  06?  J J 7,  <=  0 7 1) 

»c?  = 'S’  4 ti  - 06?) 

C 2 4 4 ; 

) 

re t urn ( i ) ; 

) 


sk i o lent ) 

I nt  ent ; 

( 

i nt  i ) 

lor  (i  - 0;  i < end  i*4) 
vr I t tlodtvioh,?!  I 

) 


C I r ( ) 

i 

i nt  i ) 

lord  = o ; i < 13?;  (44) 

1 1 f i J — • ' ; 

1 1 1 1 3 1 1 i • » ' ; 

> 


Tiov  ( t o r f roTi , n ) 
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